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Abstract 

We extensively study the exact solutions of the massless Dirac equation in 3D 
de Sitter spacetime that we published recently. Using the Newman-Penrose formal- 
ism, we find exact solutions of the equations of motion for the massless classical 
fields of spin s = i, 1, 2 and to the massive Dirac equation in 4D de Sitter metric. 
Employing these solutions, we analyze the absorption by the cosmological horizon 
and de Sitter quasinormal modes. We also comment on the results given by other 
authors. 
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1 Introduction 



The properties of de Sitter background has recently been studied [T]. In particular, it was 
found that this spacetime has several applications in different topics of theoretical physics 
such as the dS-CFT correspondence, inflation, quantum fields in curved spacetime, ther- 
modynamics of horizons [2j — [0] , since the simplicity of its metric allows us to investigate 
in detail (in some cases in exact form) many physical phenomena. Understanding of these 
phenomena in de Sitter background will be useful in more complex spacetimes, which 
have more relevance from a physical point of view. 

In order to analyze the properties of a spacetime, a method often used consists in 
studying the propagation of a classical field. In the present paper we continue to investi- 
gate the dynamics of different classical fields propagating in de Sitter metric as in Refs. 
^U]-j2I]. First, we find exact solutions of the equations of motion for several classical 
fields moving in de Sitter spacetime. In particular, for the massive Dirac equation we cal- 
culate a more simple set of solutions than those found in Ref. [21]. Using these solutions 
we study the absorption of the cosmological horizon and extend some of the results given 
in Refs. [TTH ITT] . Following these studies we compute the radial flux of a classical field 
propagating in de Sitter background when its frequency is a real number and its radial 
function is regular at r = 0. Moreover, we briefly discuss the implications of these results 
and we enumerate the differences among our results and those of other references. 

We expect that the quasinormal modes (QNMs) of the black holes will be a part 
of the signal that the gravitational wave detectors will detect in the near future |22] . 
More recently, owing to the dS-CFT and AdS-CFT correspondences [HIE], [23] • there is 
an increased interest in calculating the values of the quasinormal (QNM) frequencies in 
several asymptotically de Sitter [TEJ |21] and anti-de Sitter black holes |23l2ni- Several 
papers have studied the QNMs of de Sitter spacetime [12J — [iBJ ; nevertheless, these papers 
report different results on their existence. 1 

Motivated by this discussion in the literature on the cases in which the de Sitter 
QNM frequencies of massless fields are well defined and using the exact solutions that 
we previously found, in this work we calculate the de Sitter QNMs of a massless Klein 
Gordon field, an electromagnetic field, a gravitational perturbation and a massive Dirac 
field in 4D and of a massless Dirac field in 3D. We also comment on some of the previous 
results. 

This paper is organized as follows. In Section El we review some features of the dy- 
namics of a massless Dirac field propagating in 3D de Sitter spacetime. We improve the 
results given in Sections 3 and 5 of Ref. [T7j by using a different triad to that of the previ- 
ously mentioned reference and we find additional exact solutions of the radial equations. 
In Section |H] we obtain exact solutions of the equations of motion for the massless fields 
of spin s = 0)1)2 propagating in 4D de Sitter metric. In Section 0] we develop another 
method which allows us to find exact solutions of the massive Dirac equation in 4D de 
Sitter spacetime. Employing the new exact solutions, in each of the previous sections, we 
can calculate the absorption and QNM frequencies of some classical fields moving in de 
Sitter background. Finally, in Section we make some comments on these results. 

1 The anti-de Sitter QNMs has also been calculated in Ref. ^1] 
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In Appendix |X] we review some facts used in Section|21 In Appendix[B]we present a dif- 
ferent method for calculating the QNM frequencies given in Subsection 13.21 In Appendix 
O to compare with the solutions calculated in Sec. EJ we briefly describe the transforma- 
tions employed in Ref. [21] to find exact solutions of the massive Dirac equation in 4D de 
Sitter spacetime. 

2 Massless Dirac field in 3D de Sitter spacetime 

In Ref. we have shown that using a null triad (Eq. (3) of Ref. [T7j ) it is possible to 
find exact solutions of the massless Dirac equation when the background spacetime is 3D 
de Sitter. In the present section, we obtain exact solutions of the massless Dirac equation 
by means of a slightly different method to that used in ^7|. Employing these solutions 
we study the absorption and QNMs of a massless Dirac field when it is moving in 3D de 
Sitter metric. 

Note that in Ref. ^7] we do not consider all the cases for which it is possible to find 
exact solutions of the radial differential equations; therefore, in the following we also give 
other values for the parameters of the solutions that were not analyzed in this reference. 

In the circularly symmetric metric 



ds 2 



P{r) 2 dt 2 - Q(r) 2 dr 2 - r 2 d6 2 , 



(1) 



if we use the diagonal triad 




(2) 



and the representation of the Gamma matrices 7, 




(3) 



where <7j denote the Pauli matrices; the massless Dirac equation 



7 a V Q ^ = 7 a (e a + r a )^ = 



(4) 



becomes 



We propose a separable solution of Eq. (J""} in the form [T7J 

where m = ±|,±|,±|,... 2 Using this ansatz we find that the functions P^ and P^ 
satisfy the coupled differential equations 3 

r d iuir r/l P '\ 1 „ „ 



Qdr P 2Q \r P 

These equations are equal to Eqs. (12) of Ref. |T7]. To find exact solutions of the radial 
equations (J7J when the background spacetime is 3D de Sitter, the method employed in Ref. 
[TT] directly uncoupled the functions R\ and R2 from the system of differential equations 
J7J). We simplified to hypergeometric differential equations the second order differential 
equations obtained for R\ and R2 by means of an ansatz and a change of variable. 

To solve the coupled differential equations (J7J in 3D de Sitter metric, that is, when 

we propose that the radial functions R\ and R2 take the form (see also Ref. [2*7] ) 

R2 = 9R2, (9) 



where 



9=-T7\ ^> ( 10 ) 

22(1 — 2^ 4 



and we define the variable 2; by r = zl. 

Substituting Eqs. Q and ("IT"]) into Eqs. (J7J, we get the simplified system of differential 
equations 

where u = ul. From this system of differential equations it is possible to find that the 
function R\ satisfies the differential equation 



d 2 / 1 z \ d u 2 iuz iu m 2 



dz 2 \z P 2 dz P 4 P 4 P 2 z z 2 P 2 



R x = 0. (12) 



2 See Appendix 1X1 

3 As in Ref. |17j . for simplicity in the notation, we only write the superscript (m) of the radial functions 



when it is strictly necessary. 
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The function R 2 satisfies a similar equation. 
Now, making the change of variable ^7j, [2*U] 



1 - z 



y 



l + z 

we find that the differential equation (fT2"|) becomes 

f 



(13) 



— +( — 

dy 2 



2y 1-yJ dy 



d_ 



(J 1 — ioj 



Ay 2 



m 



22/(1 — y) y{l-yf 



Ri = 0, 



(14) 



(i?2 satisfies a similar equation). 

Making the following ansatz for Ri, 



Ri = (l- y) Bl y Cl Ri 



(15) 



with I = 1,2, and substituting into Eq. (|14j) (and into the corresponding equation for 
i? 2 ), we get that the functions Rj must be solutions of the differential equations 4 



d 2 d 
1/(1 - i/) J"2 + (°i ~ ( a i + bi + 1)2/)^— - aibi 



Ri = 0, 



(16) 



where 



a 1 = C 1 +B 1 + -- 
b 1 = C 1 + B 1 + y , 



~2~' 



a 2 = C 2 



1 ZCiJ 

5 2 + - H . 

2 2 ' 

id; 



Co + So — 



d = 2d 



1 



c 2 — 2C 2 



2' " ~ " 2 

provided the parameters B%, B 2 , C\, and C 2 take the values 

B x = ±\m\, B 2 

ft=i±!±|, c 2 



(17) 



±|m 
1 



1 iu 
- =F - ± — . 
4 4 2 



Note that in Ref. ^1 we considered that the parameters Bi and Cj take only one 
value (Eq. (26) of Ref. [T7|): however, we see in Eq. (fT%|) that these quantities can take 
two values. As a consequence, we get a larger number of solutions of the radial differential 
equation (|12j). This also happens for the massive Klein Gordon equation in de Sitter metric 
P53 and for the Dirac equation in a 3D black hole of EMDA theory [57] . 

In the rest of this section we use the variable v defined by (see Ref. |17j ) 



v 



i-y- 



(19) 



4 We notice that the symbols A, B, C, a, b, and c represent different quantities in each section of this 
paper. 
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Therefore the differential equations become 



d 2 ~ d 

v (1 - v )-^ + (c/ - {a-i + bi + 1 ) V )-^ J ~ a i b i 



Rt = 0, 



(20) 



where 



a\ — ax, 

b~i = h, 

ci = 2Bt + 1, 



02 

h 

C2 



a-2, 

25 2 + 1. 



(21) 



As is well known, the solutions of the differential equations (}2T)|) are hypergeometric 
functions [2E1- I n the following, we study in detail the case B\ = B 2 = \m\. Later we 
shall comment what happens when B\ and B 2 take other values. 

From Eqs. (|2*T|) we note that c\ and c 2 are positive integers. Thus the solutions of the 
differential equations (|2T)j) are [28J 



R 



(i) 



(2) 



+ 



2 F 1 (a I ,b I ;c I ;v), 
2 F 1 (d I ,b I ;c I ;v)ha.(v 
(ci-1 



(22) 



r(a 7 )r(& / ; 



CJ-l 

9=1 



1 



9-1 



(g-i; 



' (g/-ff-l)! 



9 =0 



Q/)g(6j) V |V(a 7 + g) + ^(6/ + g) - ^(c/ + g) - ^(1 + g) 



g!(c/) g 

-V»(5j - cj + 1) - ^(6/ - cj + 1) + + - 1 



(23) 



where = dlnr(x)/dx. 

Employing Eqs. 0, (fHjj) . and (|2*3j) it is possible to show that the radial functions that 
include i?} are divergent at r = (the origin of the de Sitter spacetime). We have no 
reasons to expect singularities of the fields at the origin of the de Sitter metric; therefore, 
in the present work we study the solutions of the equations of motion that are regular at 
this point, as was the case in Refs. jH], JD] Thus we only consider radial functions 

that include Rj. 

The parameters C\ and C 2 can take two different values for given B\ and B 2 \ then, 
it follows that it is necessary to study four cases for the possible permutations of C\ and 
C 2 , (something similar happens in the analysis made in Refs. JH], |27j ) 



Case I: 


Ci = 


Case II: 


Ci = 


Case III: 


Ci = 


Case IV: 


Ci = 



1 iu) 

2 + Y' 

1 iu 

2 + Y' 

~Y' 

iu 

~Y' 



c 2 

^2 

c 2 



1 zu) 

2 ~ Y' 

Y' 

1 zu) 

2 ~ Y' 

Y' 



(24) 
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From Eqs. (J0J), (|15p. (|22p. and using the property of the hypergeometric functions [2E] 

2 Fi(a, b; c; v) = (1 - w) c - a - b 2 Fi(c - a, c - 6; c; v), (25) 

we easily show that the functions R\ and i?2 satisfy 

\Rx\ 2 - \R 2 \ 2 = (26) 

in the four cases enumerated in Eq. (|24j). whenever Bx = B 2 = \m\ and the oscillation 
frequency co> is a real number. 

Therefore, from Eq. (37) of Ref . , the number of fermions that cross the circle of 
radius r per unit time is equal to 



dN 
~dt 



/ 2\V 2 



(m) 1 2 
2 



). (27) 



This quantity vanishes because Eq. (J26J) holds for every m when the frequency u is a 
real number. From this fact follows that the radial flux of fermions from or toward the 
cosmological horizon is equal to zero, provided we use regular radial functions at r = 0. 
When B\ = B 2 = \m\, we obtain the same conclusion in the four cases enumerated in Eq. 
()24|) . as can be explicitly verified. 

Now, we study what happens when Bi takes the value — \m\. For this value of Bi the 
parameter c x is an integer, c x < 0; thus the first solution of the differential equation (j2T?j) 
is |2H] 

rW = v 2 ^ 2 Fi(a 1 - ci + 1, h - c x + 1; 2 - c a ; u). (28) 

Taking i?^ given in Eq. ()28|) as a part of the radial function R\ y we find that it is 
regular at r = 0. The second solution of Eq. P0|) is [2H] 



^(2) = v 2|m| | 2Fl ( 5l _ ci + 1, 6i - ci + 1; 2 - ci; u) ln(v) 



+ ( — (29) 



r(oi - ci + i)r(6i - ci + 1) 



X' 



^ (5i - Cx + l) q {b x - cj + r „ ~ 

+ 2^ 1/ 9 _ - \ ^( a i - ci + 1 + g) + ^(1) - 

q=0 q ^ z Cl >i 

+ ?) - V(5i) + - ci + 1 + q) ~ i>{2 - c x + q) + ^(1 - Cx] 

Using this solution of Eq. ()20|). we can verify that the radial function Ri is divergent 
at r = 0. Thus, when B\ = — \m\, it is possible to find regular solutions at the origin 
of 3D de Sitter metric. We have an identical situation for the radial function R 2 when 
B 2 = — \m\. 

These facts imply that there are other options for the values of Bx and B 2 , also of the 
case Bx = B 2 = \m\ studied in detail in this section. The additional cases are 



(30) 



Case II: 


Bx = 


\m , 


B 2 =- 


\m\ 


Case III: 


Bx=- 


\m , 


B 2 = 


\m\ 


Case IV: 


Bx=- 


\m , 


B 2 =- 


\m\ 
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For all cases enumerated in previous equation it is necessary to consider that the pa- 
rameters C\ and C 2 can take two different values each. Nevertheless, from the conclusions 
previously obtained in this section, we believe that the same result on the absorption of 
the fermions by the cosmological horizon of the 3D de Sitter spacetime will be found in 
the cases that we do not analyze in this paper. 



2.1 Quasinormal modes 

The QNMs of an asymptotically flat black hole are defined as the set of oscillations 
of a classical field that satisfy, at the horizon and far from it, the following boundary 
conditions: (i) the field is purely ingoing at the horizon of the black hole; (ii) the field 
is purely outgoing at infinity [22] • If the black hole is not asymptotically flat, then it is 
possible to impose other boundary conditions far from the horizon; see for example Refs. 
|26j . |27| . [2H|, [20] • For static and stationary black holes the QNMs for different fields 
have been studied mainly using numerical methods |22j . 

An interesting problem is to investigate the analogous of the QNMs of black holes 
for the fields that are propagating in de Sitter metric. Clearly, the definition of the de 
Sitter QNMs must be different from that used in the case of a black hole. In this paper, 
following Ref. [12] (see also Refs. [T3], [Ej), the de Sitter QNMs boundary conditions are: 
(i) the field is regular at the origin; (ii) the field is purely outgoing near the cosmological 
horizon. 

Using the exact solutions of the Dirac equation previously found, we can calculate the 
QNM frequencies of a massless Dirac field propagating in 3D de Sitter background. In 
the following we study in detail the case B x = B 2 = \m\, G\ = \ + y> an d C 2 = y- 

By the property of the hypergeometric functions 



, , x T(c)T(c — a — b) , . 

2 F 1 {a, b- c; v) = — ^- 2 F X a, 6; a + b + 1 - c; 1 - v) 

T(c — a)T(c — b) 

+ r(c) 5/°xl?r C) (l - vy-^F^c - a ,c-b-c+l-a-b;l-v), (31) 

r(o)r(6) 

which is valid when c — a — b is not an integer, we can write the radial functions R\ and 
R 2 that are regular at the origin in the form 

Rl = ^h-i ( (1 _ ^mnci-ai-h) (32) 

V2 \ v r(c! - di)r(ck - h) 

x 2 FM X , b x] 1 + di + h - ci; l-v) + (l- v)l- c ^ {c ~ l)T{dl + bl ~ C ~ l} 



r(d 1 )r(6 1 ) 

x 2 Fi(ci - di, ci - b\; 1 - di - b\ + ci; 1 - u)| , 



and a similar expression for R 2 which can be easily found from Eq. (j32j) substituting the 
subscript "1" of the parameters by "2". 

According to the definitions previously given, the following relation holds 

l-v = exp(-m, (33) 
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where 

stands for the usual tortoise coordinate which satisfies r* — > as r — > and r* — ► oo as 
r -> Z [Hj. 

Employing Eqs. © and f)33|) . we get that the second term in curly braces in the right- 
hand side of Eq. (|32j) represents an outgoing wave, whereas the first term in curly braces 
represents an ingoing wave. Thus, in order to satisfy the definition of the QNMs in de 
Sitter spacetime, we must cancel the first term in curly braces of Eq. One way to 

achieve this cancellation is to find the poles of the terms T(c x — d x ) and T(c\ — b\) [12 , 
[T3] . These poles are localized at [2E| 

c\ — di = —ni, and c\ — bi — —ni, (35) 

respectively, where n\ = 0, 1, 2, . . . 

From the values of di, b\, and c\ given in Eq. ()21j). the previous equations yield 

\m\ = —rii, \ m \ + - — i£o = —fii, (36) 

respectively. We cannot satisfy the first relation in Eqs. (}3T)j) (that is, T(ci — bi) has no 
poles), but the second relation in (jSTj) implies 

1 

iuo = \m\ + — + rii. (37) 

Note that the term r(ci) does not have poles, whereas the poles of the term T(ci — 
di — bi) are localized at 

ci — di — b\ — —n 1 . (38) 
Using the values of Si, b\, and ci we find that the previous equation becomes 

1 . , 

iu = rii — -. (39) 

Since the quantity m is a half- integer (see Appendix EJ, the expression in the right- 
hand side of Eq. ()37j1 is an integer, whereas the corresponding term in Eq. ()39)1 is a 
half-integer; therefore they cannot be equal, and there is no pole cancellation in the first 
factor of the expansion of the radial function (|3~2j). (see Refs. [T2], [H], [El)- We notice 
that for the frequencies iu given in Eq. (|3"7j) the quantity ci — di — bi is not an integer, 
that is, if the quantity m is a half-integer, then there is no pole cancellation. 

By analogy with the previous case, the radial function R 2 will be purely outgoing as 
r — > Z if either of the conditions 

c~ 2 - S" 2 = -ni, or c~ 2 -b 2 = -n x , (40) 

holds. In previous equations we can satisfy the first condition but not the second one. 
The QNM frequencies obtained are equal to those given in Eq. (pT7j). If Bi = B 2 = \m\, 
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then for the other possible values of C\ and Ci (see Eq. (J24|0 . we also find the QNM 
frequencies ("3*7)) . 

Taking into account that the exact solutions of the massless Dirac equation studied in 
the present section have a harmonic time dependence of the form e~ lujt , we find that for 
the QNM frequencies (13 7|) the amplitude of the field diminishes with time. 

When a massless Dirac field has some of the frequencies (j3*7|) . an explicit calculation 
of the radial flux (J27)) shows that it is different from zero. Therefore the modes with 
frequencies (f3~7)) are QNMs. 

The frequencies ([*"?]) are equal to those given in Eqs. (39) and (40) of Ref. [T*""j when 
the mass of the Dirac field is equal to zero and if the quantity £, which appears in Eqs. 
(39) and (40) of Ref. [TH], i s a half-integer greater than zero (Eq. (39)) or smaller than 
zero (Eq. (40)). The exact values allowed for I are not explicitly found in Ref. [TKj . 

The existence of the QNM frequencies for a massless Dirac field implies that for this 
case the conjecture of Ref. [H] holds. For a massless classical field, this conjecture proposes 
that there are well defined de Sitter QNMs only in odd spacetime dimensions. (See below.) 

As previously discussed in Ref. [14] . the boundary conditions imposed in Ref. [TTJ on 
a classical field propagating in de Sitter spacetime to calculate their QNM frequencies are 
not equal to those imposed in Refs. [12], [13], |Ej [E] (and in the present paper). In 
[TT] . the de Sitter QNMs boundary conditions are: (i) the field has zero flux at r = 0; (ii) 
the field has a outgoing flux different from zero near r = I. Moreover, in Ref. |JT| only 
real frequencies are studied, therefore it is deduced that the de Sitter QNM frequencies 
do not exist. As is well known, the QNM frequencies are complex numbers. For these 
reasons we do not compare our results with those of Ref. [TT] . 



The behavior of the classical fields in AD curved spacetimes can be analyzed with an 
efficient method that employs the Newman- Penrose formalism [3*1] . [3*2] . [3*3*] in order to 
get the equations of motion for the fields in a suitable form. In this section we study 
some aspects of the dynamics of massless fields propagating in de Sitter metric in four 
dimensions using this method. 

Henceforth we use the null tetrad 5 [3*IJ, [32] 



3 Massless fields in AD de Sitter spacetime 




(41) 



V 



1 



(0, 0,1,2 csc 9) , 



m' 



V2r 



1 



(0,0,1,-i csc#), 



m 



V2r 



5 In this work z* stands for the complex conjugate of z. 
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where P 2 was previously defined in Eq. (jSJ) . The tetrad (jH)) is different from that used in 
Refs. [211], jH| : however, we employ (|HJ) because the necessary algebraic steps are simpler 
than in the tetrad used in j20j, [2T| . 

The spin coefficients for the tetrad ([4*T]) are equal to 6 

K = a = X = U = T = 7f = 0, 



p = H 

e = 7 : 



P 



a 



r 

cot6> 
2V2r 



(42) 



For the Weyl tensor, the traceless part of the Ricci tensor, and the scalar curvature we 
get 



M> = q? 1 =\& 2 = v]>3 = v]/ 4 = 0, 

$00 = $01 = $02 = $11 = $12 = $22 = 0, 

f 



(43) 



As is well known, in a type D or simpler spacetime (in particular de Sitter back- 
ground), the equations of motion for the massless fields of spin s = 5, 1, 2 admit separable 
solutions [32], [33]. The equations for the field components invariant under gauge 
transformations and infinitesimal rotations of the tetrad are given by (221, E31 



{[D - (2s - l)e + e* - 2sp -p*][A + ji- 2s 7 ] - (2s - l)(a - 1)* 2 
-[5 + r - 2sr - a* - (2s - 1)0\[S* +tt- 2sa]}xt = 0, 

for s = -l-i, +1, +2, and 

{[A - (2s + 1)7 - 7* - 2sp + /}*][£> - 2se - p] - (2s + l)(s + 1)* 2 
-[5* -t* - 2sfr + (3*- (2s + l)a}[5 - r - 2sfi]} X ; = 0, 



(44) 



for s = -|, — 1, 



-2. 



In AD de Sitter metric, using the tetrad (J41j) and proposing the ansatz 



-iujt 



Rt(r) s Y< 



(45) 



(46) 



where the symbols s Yjm stand for the spin- weighted spherical harmonics [35] , we can show 
that Eq. (J44j) simplifies to 



dz 2 



2z 



2(5 + r 



1 



d 

dz 



Cj 2 



+ 



2iujs 



2s 2 + s + 2 



s 2 z 2 



[1-z 2 ) 2 z(l-z 2 ) 2 
(j + s)(j-s + l) 



+ 



2s 



1 



z^ 



(1 



z 



2\2 



Z 2 (l 



Rt = 0, 



(47) 



6 For the spin coefficients /1 and tt in the original Newman-Penrose notation in the present paper 
we use the symbols jl and tt respectively. 



11 



where u and z are defined as in the previous section. 

Employing the variable y defined in Eq. ()13j). the differential equation (J47)) becomes 



{dy 2 + 



2(s + l) 2(s + l) 



1-2/ 



1 + 1/ 



d 



£ + is) 



j(j + l)-a(s + l) jO' + 1) - s(s + 1) 



1/(1-2/) 



2/(1-2/) 
(s + 2)(s+ 1) 



(1-2/) 2 
(s + 2)(s + l) 



y(l + y) (1 + 2/) 2 

Proposing that the function Rf takes the form 7 



0. 



(l+2/) A+ (l 



2/)*V + H 



with the parameters 



A + 
B + 

C + 



-0' + 1*1 

M 

2 ' 



— 

2 



2 ' 2 ' 



(4f 



(49) 



(50) 



we find that the function i?t must be a solution of 



b + + l)v)— -a + tf 
av 



with t> 



y, as in the previous section, and 



, „, _ , iCj 3\s\ 

a+ = C + + B + + 1 + — + -Li, 



c + = 2(fi+ + \s\ + 1). 



T + 2 ' 



(51) 



(52) 



Taking into account the possible values for B + given in Eq. ([50)1 . we see that c + is 
a positive (for Z? + = j — \s\) or negative integer (for B + — — (j + \s\ + 1)). As in the 
previous section, in both cases, it is possible to find solutions of Eq. (J5lf that lead to 
regular radial functions at the origin of de Sitter spacetime. Moreover, in both cases 
there are solutions of Eq. (|5*T]) that lead to radial functions Rf which are divergent at 
r = 0. 



7 In Ref. [23j was shown that the radial differential equations of the massless fields of spin s — 0, h, 1,2 
moving in de Sitter background can be simplified to differential equations of hypergeometric type. The 
tetrad used in Ref. [201 is different from that used in the present paper, and we believe that the method 
used here to make the algebraic simplification is even simpler, (see also Section 0}. 
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In the following we only study in detail the case B + = j — \s\. For the other case (when 
B + — — (j + \s\ + 1)) we expect to find identical physical results. If B + — j — \s\, then 
the quantity c + is a positive integer. Therefore the solutions of the differential equation 
(EEJ) are J2H1 



Rt {1) = 2 F 1 (a + ,b + ;c + ;v), 
R+W = 2 F 1 (a + ,b + ;c + ;v)\n(v) 



(53) 



+ 



T(a+)T( 



q=l 



a + + q) + 4){b + + q)- ^(c + + q) - V>(1 + ?) 



9=0 

-V^(a + - c+ + 1) - ^(6+ - c + + 1) + V(l) + ^(c + - 1 



(54) 



The second solution Rf^ of Eq. (j51j) leads to a radial function which is divergent 
at the origin. Therefore in the rest of the present section we only consider the first solution 
R,t of the differential equation (|51|). which hereafter we simply denote as Rf. 

A similar procedure with 



-iuit 



R7(r)„Y< 



s ± ]my J i f / 1 



(l + y) A (l-y)» y<< R 
1 + ls 



X s 
R' s 
A' 

B' 
C 



simplifies Eq. ()45|) to a differential equation of hypergeometric type for the function R 

d 2 



(55) 




v(l - v)— + (c~ - (aT + 6" + l)u) 
dtr 



d 

df 



a b 



R: = o, 



(56) 



where 



G- 



B- 



B~ 



IUJ \s\ 

Y + Y 



lid 

Y 



3|s| 



(57) 



c~ = 2{B- + \s\ + 1). 



Depending on the value chosen for B~, the parameter c~ is a positive or negative 
integer. In both cases there are solutions which are regular at the origin of de Sitter 
metric. In the rest of the present section, we only study in detail the case B~ = j — \s\, 
that is, when the quantity c~ is a positive integer. 
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If B = j — \s\, then making a similar analysis for R s , we find that the solution of 
the differential equation (J5HJ) 

R- = 2 F l (a" 1 b"]c-]v) (58) 

leads to a radial function R~ which is regular at r = 0. 

As in Section |21 for given B + and B~ , it is necessary to note that the parameters C + 
and C~ can take two values each. This fact implies that we have the cases (see also Refs. 



Case I: 


C + 


\s\ 

~ ~Y 


iuo 

+ r 


c- 


= + T 


iuo 

+ ¥• 


Case II: 


c + 


\s\ 

~ ~Y 


iuo 

+ y 


c~ 


\s\ 

~ ~Y 


iuo 

' T' 


Case III: 


c + 


= + i 


iuo 


c- 


= + T 


iuo 

+ ¥• 


Case IV: 


c + 


\s\ 

= + i 


IUJ 

~ y 


c~ 


\s\ 

~ 2 


iuo 

~ Y' 



(59) 



If the frequency uo is a real number and the radial functions are regular at the origin 
of de Sitter spacetime, when B + = B~ — j — \s\ it is possible to show that the relation 

\Rf\ 2 - \r;\ 2 = (60) 

holds in the four cases enumerated in Eq. for s = |, 1, 2. 



3.1 Flux calculation 

In this subsection we calculate the radial flux of the massless fields of spin s = |, 1, with 
regular radial functions at the origin of de Sitter metric. The analysis for a scalar field 
(s = 0) was published in Refs. [TU], [TTj . 

If the spin of the massless field is s = |, then we use that the radial component of the 
conserved current for a massless Dirac field in the tetrad (}4~T|) is equal to [121 

^ = ^(W-ixirV (6i) 

V 2 2 / 

Therefore the number of fermions that cross the sphere of radius r is given by 
ON r:: 



dt 



pZTT PIT / \ 

/ / J r ^jd^d0ocPr 2 (| J Rl| 2 -| J R 1 | 2 . (62) 

JO JO V 2 2 / 



This quantity is equal to zero due to Eq. (|5U|) . 

Using Eq. (339) given in Chapter 5 of Ref. [321, we find that the component T\ of 
the energy-momentum tensor of an electromagnetic field propagating in 4D de Sitter 
background is equal to 

T r t =^(\xl\ 2 -\xtn (63) 
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when we use the tetrad (jUJ). 

Integrating Eq. (J63)) over the coordinates 9 and <p, we find that the energy flux of an 
electromagnetic field takes the form 

^ocr 2 P(|^| 2 -|i? r | 2 ), (64) 

which is equal to zero owing to Eq. (|6U|) . Thus for a massless Dirac field and an electro- 
magnetic field moving in de Sitter spacetime Eqs. (|62j) and (|64|) show that their radial 
flux is equal to zero whenever the radial functions are regular at r = and the frequency 
u is a real number. This conclusion holds in the four cases enumerated in Eq. (joTJj) . 

However, we have not been able to prove the equivalent result for the gravitational 
perturbations. By using the Isaacson stress-energy tensor j^S], we found a complicated 
expression for the energy momentum flux at any point in de Sitter metric. Although we 
can use Eqs. (jH|) and (J4l)|) to describe the dynamics of a massless field of spin s = | 
(Rarita-Schwinger) in 4D de Sitter spacetime jSZj, we do not know an expression for its 
energy-momentum tensor that allows us to calculate its radial flux. 



3.2 Quasinormal modes 

By the definition of de Sitter QNMs previously given in Subsection 12.11 it is clear that the 
solutions of the equations of motion for the massless fields found in the present section 
allow us to calculate these by means of a similar procedure to that of Refs. [T2] . [T3j . jHj. 
In the following we study in detail the case B + — j — \s\ and C + = y — -y. F° r these 
values of the parameters B + and C + , the regular solution at r = is 

, . . . iu) \s\ 

R+ = (2 - t>) 1+ IV-l s l(l - f)T-T 2 F 1 (a + , b + ; c+; v). (65) 

Employing the property of the hypergeometric functions ()31|) valid when c — a — b is 
not an integer [28J, we can write the function i?+ of the previous equation in the form 

JH 2-,)'^'4(i-.if-? r(tt)r|tt - aMt) 



T(c+ - a + )T(c^ 



iu \s\ 

x 2 Fi(a + , b + ; a + + b + + 1 - c + ; 1 - v) + (1 - «)" T+! (66) 

r(c+)r(a+ + b+ - c+) , , , , , 

x v — r— y— -r ^ 2 ^i(c + - a + , c+ - b + ] c+ + 1 - a + - b + ; 1 - v) 



From Eq. (|46|). we note that the massless fields studied here also have an harmonic 
dependence on time e~ luJt ; therefore the first term in curly braces of Eq. (|66[) represents 
an ingoing wave as r —>■ I, whereas the second term represents an outgoing wave as r — >• /. 
To satisfy the boundary conditions of the de Sitter QNMs, we look for the poles of the 
terms T(c + — a + ) and T(c + — b + ) (see Subsection 12. 1|) . These poles are localized at |2E] 

c + — a + = —rii, c + — b + = —n>i, ni = 0, 1, 2, . . . (67) 

Hence the quantities c + — a + and c + — b + must be negative integers, but from the 
values of a + , b + , and c + given in Eqs. (|52*|) . we see that if c + — a + — b + = —{iu — \s\) is not 
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an integer, then the quantities c + — a + = j + 1 — \s\ — (iu — \s\) and c + — b + = j + 1 + \s\ 
cannot be negative integers. That is, we cannot satisfy the conditions (JB7J). Therefore, if 
the quantity c + — a + — b + is not an integer, then there are no QNM frequencies for the 
massless fields of spin s = |, 1, 2. 

In Refs. [12], ^3] this fact was not noticed; moreover, in these papers there is the 
so called pole cancellation, which implies that a massless minimally coupled scalar field 
has no well defined QNM frequencies when it is propagating in 4D de Sitter spacetime. 
Something similar happens in Ref. for the gravitational perturbations if the dimension 
D of de Sitter spacetime is even, D > 4. (See Appendix iBl) 

We stress that if the factor T(c + — a + — b + ) has poles, then Eq. (|ST|) is not valid 
When c + — a + — b + = —n, n = 1, 2, . . . , instead of Eq. (Ell) , we must use 8 j2Hj 



n , i \ T(a + b - n)T(n) ^-r, 1 ^ i a ~ n )ci{b ~ n) q . 



_L?ffl!^ E J^ (1 _ rt . Wl _„ ) 

1 (a — njl (o — n) + q)\ 

- ip{q + 1) - tp(q + n + 1) + V(a + ?) + ^(6 + ?)]■ (6? 
From Eq. (J68j) . we find that the function i?+ takes the form 

R+ = (2- v) 1+l % 3 - ls \l -v)2 

r(q+ + b + - n)Y{n) \-n\^ ( Q+ ~ n )i( b+ ~ n )i n ^ 



r(a+-n)r(6+-n) ^ g!(n + g)! v Jl v ; rvy y 
- ^(g + n + 1) + ip(a + + q)+ ip(b + + ?)]}■ (69) 

In this expression for Rf, the second term in curly braces represents an ingoing wave as 
r — ► I; then, in order to have a purely outgoing field as r — > /, it is necessary to satisfy 
the condition 

a + — n = —rii, or b + — n = — m. (70) 

From the values for a + , 6 + , and c + given in Eq. (J52|) we cannot satisfy the first condition 
in Eq. (J7U|) . but from the second condition we find that the QNM frequencies are equal 
to 

iu = j + l + m. (71) 

Making a similar analysis for R~ , when B~ — j — \s\ and C~ = y + y , we a l so find the 
de Sitter QNM frequencies (J7TJ). We point out that for the frequencies (J7TJ) the amplitude 
of the field decreases with time. For massless fields of spin s = ~, 1 with frequencies given 

8 If c + — a + — b + = n or c + — a + — b + = 0, n = 1, 2, 3, . . . , then there are no QNM frequencies because 
the amplitude of the purely outgoing fields increases with time or the fields have no time dependence. 
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by Eq. ([71)1 . an explicit calculation of the radial flux shows that it is different from zero 
for each mode. 

The dynamics of the gravitational perturbations propagating in D-dimensional de 
Sitter metric (D > 4) has been studied in Ref. using a different set of exact solutions. 
Moreover, the authors of Ref. [T3. conjectured that the massless classical fields moving 
in de Sitter background do not have QNM frequencies when the spacetime dimension is 
even. The results presented here show that this conjecture is not valid in four dimensions 
for the massless fields of spin |, 1, 2. In particular our result for the QNM frequencies of 
the gravitational perturbations propagating in 4D de Sitter metric is different from that 
of dU. 

To finish the present section we note that it is possible to consider other permutations 
for the values of B + and B~ , different from the combination studied in detail here; for 
example, the combination B + — j — \s\ and B~ — — (j + \ s\ + 1) (and other two options). 
Nevertheless, we believe that in the other cases the same results hold for the value of 
the cross section of the cosmological horizon and for the values of the de Sitter QNM 
frequencies. 

4 Massive Dirac field in AD de Sitter spacetime 

When the background spacetime is 4D de Sitter in static coordinates, some time ago 
Otchik found exact solutions of the differential equations system for the radial functions 
of the massive Dirac equation [21] (see Appendix EJ) • Here we present a different method 
to get exact solutions of the radial equations employing hypergeometric functions. 

For the massive Dirac equation we make the separation of variables in the tetrad (|4*T|) . 
which is different from that used in Ref. [21]; nevertheless we obtain the same system of 
differential equations for the radial functions. 

We write the massive Dirac equation using the Newman- Penrose formalism in the form 

IS], 132] 

(D + e- p)F 1 + (6* + 7f - a)F 2 = i^G x , 

(A + p. - 7 )F 2 + (5 + (3 - r)Fi = i^G 2 , (72) 
(D + e* - p*)G 2 -(8 + r- a*)G x = i^F 2 , 
(A + PT- 7 *)Gi - (5* + r*)G 2 = i^F u 

where M = a/2/x* is the particle mass. 

Employing the tetrad (j4*T|) in 4D de Sitter metric and making the ansatz 

Fi = e-^_ L Y jm (e, 0)A(r), F 2 = e'^iY^, <P)F 2 (r), 

G x = e- <wt _il} m (0, 0)Gi(r), G 2 = iY jm {6 , <f>)G 2 {r) , (73) 

where F 2 = G% and Fi = G 2 , it is possible to show that the massive Dirac equation in 
4D de Sitter background simplifies to the system of differential equations for the radial 
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functions [2T], EH 



P-r - tt - 777^7 + - )Fi + — iM F 2 



dr P 2P1 2 r J \ r 

d : iuo r P\ ^ ( j + \ 



/> d; + p-2PF + 7j^HV +tM J F ' = °- (74 > 

where j is a half-integer, j > ~. 

Making the changes F\ —* R\, F% —R 2 , r — zl, u — ul, and M = M/ we find that 
the coupled differential equations (J73j) become 

1 - ^ — — — -y ^ +^^- 1 «x 

dz (1-32)* 2(1-Z 2 )2 Z ' 



i + 1 



z 



-iM) R 2 



+ z -^+^^)r, 

dz (l-Z 2 )* 2(l-2 2 )5 



z 



iM R x . (75) 



Proposing the ansatz 



where 



z(l -z 2 ) 

we can simplify the system of equations (J75J) to 



i?i = hR\, 

R 2 = hR 2 , (76) 



h= - 1 -t , (77) 



.2a d ^ \ 6 _ M + 2 



'l-z'Y--- -r \R 1= r-^-iM R 2 , 

dz (1 _ z 2 )i J \ z ) 

d^ (l-z 2 )sj y 

The system of differential equations (|7Hj) is equivalent 9 to that given in Eq. (71) of Ref. 
[TT] (from this system we can easily find Eqs. (4) of Ref. [2T], as discussed in Ref. |T7|). 
In Appendix O we briefly describe the method used in Ref. [2TJ to find exact solutions 
of the system of differential equations (f?8*J) employing hypergeometric functions. In the 
following, we explain a different method to obtain exact solutions of the coupled equations 
(1781) . which we believe is simpler than the one used in 

identifying j + h with m and M with jl. 
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First, we make the ansatz (compare with that given in Eqs. (5) and (6) of Ref. |2*T] . 
which we write in Eqs. (|118|) and (|119j) of Appendix 

fli = (l-^)-i(l -/ 2 ), 

R 2 = (l-z 2 )-^l + z) 1 Hf 1 + f 2 ). (79) 

Substituting these expressions for the functions R\ and R 2 into the differential equations 
dZHJ) yield 

(i-^ 2 )^(a+/ 2 )+Q+^) (/i+/ 2 ) 

- - z){h - f 2 ) - iM{\ - z){h - h) = 0, 



(l- z ^—( fl -f 2 )- I +& (/i-/ 2 ) 



dz 



(80) 



(1 + + / 2 ) + zM(l + + / 2 ) = 0. 



From previous equations, we find that the functions fi and f 2 must be solutions of 
the coupled differential equations 



(I-, 2 ) 



dz 



l^A-iMz 



- iMz 



fi = - 
h = - 



- + uj-{j + -)+iM 

- + iu+{j + -)-iM 



/a, 
A- 



(81) 



We easily find from this system of differential equations that the functions f\ and f 2 
satisfy 



(1-z 



2\2 



dz 2 



- 2z(l - z 2 ) — + {l-z 



dz 



dz 2 



dz 



3 + 



+ iM 



(^i-iMz^ Ja = o, 



(1 - Z 2 Y-^r - 2z(l ~ Z 2 )-^ - (1 - Z 2 ) 



+ iM 



3 + h 



-iMz) }f 2 = 0. 



(82) 



Making in Eqs. (|82jl the change of variable x = z 2 , and taking fi and f 2 equal to 

f 2 = x c %l-x) B >f 2 , ( 



53) 
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where 



Ci 



B 1 = ± 



2 ' 
1 fl 



C 2 



5 ± 0' + 1) 



2 V2 



iuj 



B 2 = ± 



1 / 1 



2 \2 



+ ia> 



we find that the functions f\ and f 2 are solutions of the differential equations 

A = o, 



d 2 d 
x(l - x )-^ + i c i - ( a i + bi + _ a ibi 



where / = 1, 2, and the parameters ai,bi, and cj are equal to 



r, ^ 1 iM 



D 1 iM 
a 2 = ±s 2 + ^2 + —, 



b 1 = B 1 + Ci 
ci = 2d + -, 



iM 



B 2 + C 2 + 



iM 



c 2 = 2C 2 



54) 



(86) 



Therefore the functions /i and / 2 are solutions of hypergeometric type differential equa- 
tions. 10 

As in previous sections, we have several sets of possible values for the parameters of 
the exact solutions depending on the values chosen for B\, B 2 , C\, and C 2 . (See previous 
sections and jT^], |2Tj.) 



4.1 Flux calculation 

The exact solutions of the massive Dirac equation previously found allow us to study the 
absorption of a massive Dirac field by the cosmological horizon. As is well known, 11 Eq. 
(16 1|) also determines the radial component of the conserved current for a massive Dirac 
field propagating in 4D de Sitter spacetime [32] • Therefore, from Eq. ([B2~)l. we must study 
the quantity |i?i| 2 — |.R 2 | 2 to calculate the radial flux of a massive Dirac field. Henceforth 
we analyze in detail the case B\ = B 2 = |(| + iu), C\ = |(j + |), and C 2 = C% + \. 
Employing Eqs. (|75jl . (J79j) . and (J83j) we can show that the functions R\ and R 2 satisfy 

o~4Ci-l 

\Ri\ 2 - \R2? = r {(/i + / 2 )(/r + ft) - (1 - x)f 2 ft}. (87) 

(l-^ 2 )2 

Thus it is necessary to analyze the terms in curly braces of Eq. (jHZj) . For the case studied 
here, it is easy to show that the radial solution that is regular at the origin of the de Sitter 
metric is the function that includes 

/i = 2 Fi(ai,&i;ci;x), (88) 

10 We point out that there are many papers in which the properties of the exact solutions of the massive 
Dirac equation in 4D de Sitter spacetime are investigated, see for example Refs. EU, EHh [HSj- 
11 Making appropriate changes in the notation. 
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and for this case, from Eqs. (|%Tj). we get that the function f 2 is equal to 

h = 2F 1 {a 2 ,b 2 ;c 2 ;x) = 2-^1(01, h + l;ci + l;x), (89) 

Ci Ci 

which is regular at r = 0. 

Equations (JBEI) and QSfy yield 

b* 

fi + h = —2^1(01, 61; ci + 1, x), 

Cl 

A/ 2 * = ^ir(l - *)~ ViK 6 i; ci + 1, x)| 2 , (90) 
c i 

hence 

(A + A)(A + ti) - a - *)/ 2 £ = 0, (91) 

for all r G (0,/). 

Taking into account Eqs. (|62j) . (|87j) . (|9Uj) . and (|91j) we get that the radial flux of the 
massive Dirac field from or toward the cosmo logical horizon is zero, if we use the regular 
radial solutions at r = and the oscillation frequency u is a real number. This result 
is similar to that found in the previous sections for the massless fields of spin s = |, 1 
moving in 4D de Sitter and for a massless Dirac field propagating in 3D de Sitter, (see 
also Refs. [TO], [UJ, HE). 

From the result that the cross section of the cosmological horizon is zero, in Ref. 
[TU] was asserted that it does not emit Hawking radiation (see Introduction of 10J), an 
affirmation that contradicts the well established thermodynamical properties of de Sitter 
spacetime. We note that in Refs. ^UJ, (201, |2J| a value different from zero was found 
for the absorption probability of de Sitter horizon, but in these references only radial 
functions that are divergent at r = were studied. We believe that the implications of 
the results of Refs. [TU], [TT] and this paper on the emission properties of the de Sitter 
horizon must be investigated further. 

An alternative is that our result follows from the assumption of thermal equilibrium 
made throughout this work and in Refs. [TU], [TT], [T7j, (in our analysis the de Sitter 
cosmological horizon does not expands or shrinks). 



4.2 Quasinormal modes 



We can calculate the de Sitter QNMs of a massive Dirac field employing the exact solutions 
found in this section. When B\ = -(^ + ioS) and G\ = \{j + \) from the property of the 
hypergeometric function given in Eq. (|31|) , we can write the radial function R\ of Eq. (|76p 
in the form 



Ri 



h (1 - z)2Z j 



Cl 



(92) 



'l-z 2 



X 



r(ci + 1) 
r(ci 



ir( Cl 

+X2 



6l) 
- dl 



r(d - ai - h) 

Wei 



r(ci + 1 — ai) 



2 F 1 (ai,b 1 



2 Fi(ai, 61; ai + 61 + 1 - ci, 1 - x) 



1; a x + bi + 1 - Ci; 1 



x) 



x)^~ 
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r(ci + i)r( ai + 6i-ci) r 7-1 \ 

H T-i / ^77 , 1N 2-Pi (ci - a u ci - 61; c x + 1 - a x - 6 a ; 1 - x) 

r(ai)r(6i + 1) 



1 

+x2 2 Fi(ci - ai + 1, ci — 61; ci + 1 — ai — 6 X ; 1 — x) 



In Eq. ()92|) . as r — >■ I, the term with the factor (1 — x) 2 ^ 2 ) represents an ingoing 

wave, whereas the term with the factor (1 — x) 2 ^ 2 > represents an outgoing wave. 
Therefore, in order to have purely outgoing oscillations as r — > I, it is necessary to satisfy 
the condition 

ci — ai = — ni, or ci — bi — — n>i, % = 0, 1, 2, . . . , (93) 

which yield 

iu = j + 2ni - iM, iu = j + 2n x + 1 + iM, (94) 

respectively. 

We easy see that the factor T(ci + 1) does not have poles, while the factor T(ci — at — 61) 
has poles at 

iu = n\ — -. (95) 

From Eqs. (|9~%j) and ()95j) we get that these values for iu cannot be equal if M 7^ 0, hence 
there is no pole cancellation. We note that for the frequencies (J94[) the quantity c\ — a\ — b\ 
is not an integer. 

We also find the values ()94j) for iu when we analyze the radial function R 2 by means 
of a similar method (B 2 = \{\ + iu) and C 2 = |(j + |) + | correspond to the case studied 
in detail). 

An explicit calculation of the radial flux ()62j) shows that it is different from zero if the 
frequency of a massive Dirac field is equal to some of the frequencies f!94j) . and we note 
that for these frequencies the amplitude of the field decreases with time. Therefore the 
oscillations with frequencies equal to those given in Eq. (}9"4*|) are QNMs of a massive Dirac 
field moving in 4D de Sitter spacetime. 

If M = 0, then the QNM frequencies (jUj) are equal to those of Eq. (|7T|) for s — |. 
Moreover, the frequencies ()94j) coincide with those given in Eqs. (50) and (51) of Rcf. [15 , 
if the quantities k + and k_ are positive and negatives integers, respectively. In Ref. [T3] 
the method used to calculate the QNM frequencies is different (and more complicate) to 
that applied in this section. 

An analysis of the QNM frequencies (|9*lj) show that, at least in four dimensions, the 
mass of a Dirac field does not have an inferior limit so that the QNM frequencies are well 
defined. As shown in Ref. [TS] (see also JH]), for a minimally coupled scalar field there is 
an inferior limit for its mass in order to have QNMs. 

We stress that when the frequency of a Dirac field is equal to iu = j — iM (iu = 
j + 2n\ — iM when n\ = 0), the method previously used to find the exact solutions of the 
massive Dirac equation is not applicable. Nevertheless, a detailed analysis of this case 
shows that the solutions found above hold replacing iu by j — iM in the values of the 
different parameters of the solutions previously given. Moreover, if the frequency of a 
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Dirac field is iu = j — iM, an explicit calculation of the radial flux (|62j) prove that it is 
different from zero. 

In the present section we find exact solutions of the Dirac equation in 4D de Sitter 
background, but taking into account the results of Ref. [T7j we deduce that the exact 
solutions of the system of differential equations (|74|) obtained in this section lead to exact 
solutions of the massive Dirac equation in 3D de Sitter metric, provided we make the 
identifications mentioned in footnote El Thus, making these identifications, the frequen- 
cies (194)) are also QNM frequencies of a massive Dirac field propagating in 3D de Sitter 
spacetime. These values for the QNM frequencies are equal to those given in Eqs. (39) 
and (40) of Ref. jT5] if the quantity £ is a half-integer. To calculate these frequencies the 
authors of Ref. ^H] use a more complex method. 

The method exploited in the present section also produces exact solutions of the 
massless Dirac equation when we take the limit /x — > 0. A quick look shows that the 
solutions found in the present section are different from those studied in Section El of the 
present paper. We have not been able to find the relation between the two sets of exact 
solutions. 

We believe that the usefulness of the exact solutions previously found to study other 
physical phenomena occurring in de Sitter metric must be investigated further. 

5 Discussion 

In the present paper we studied the absorption and QNMs of some classical fields propa- 
gating in de Sitter spacetime. If the frequency u is a real number and the fields have no 
divergences at the origin, the result found here and in Refs. [TU] . jTJ], ^7] for the radial 
flux of a classical field moving in de Sitter background is the same (zero). We obtain this 
value for the radial flux studying the behavior of different fields, and we believe that the 
implications of this result must be analyzed further. An analogous fact is that the event 
horizon of an extreme BTZ black hole exhibits a vanishing absorption cross section, as 
shown in Refs. |4*U] . 

Our results and those of Refs. IIS] suggest that for a massive classical field moving 
in de Sitter metric there are always QNMs. According to Ref. |14j . if the mass of the field is 
zero, then de Sitter QNMs exist only when the spacetime dimension is odd, but in Section 
13.21 fsee also Appendix [BJ), we study an example with well defined QNM frequencies when 
the spacetime dimension is even and the fields are massless. Thus the affirmation of Ref. 
[T4"] on the existence of the de Sitter QNMs is not generally valid. Is it possible to get 
other massless fields with well defined QNM frequencies in de Sitter metric? (different 
from those found in Refs. [Hj, [T5], and the present paper). We believe that this point 
must be studied further. 

For massless fields the QNM frequencies found in this paper are purely imaginary as 
those found in Refs. [Hj, j2JJ, EH], EU, 02], and we notice that for the QNM frequencies 
calculated in the present paper the amplitude of the field diminishes with time. 

We found that several parameters of the radial functions can take two values; as a 
consequence there are many possible combinations for the parameters depending on the 
values chosen. Nevertheless, some of our results indicate that the values chosen for the 



23 



parameters do not change the physical conclusions obtained. (See also Refs. [E1I2ZI-) 

An interesting question is to study the relevance of the results reported in the present 
paper for the dS-CFT correspondence [5], [E], [21] • A natural extension of the present 
paper is to investigate if its results hold for the D- dimensional de Sitter spacetime. Finally, 
is there any relation between the purely imaginary de Sitter QNM frequencies and the 
algebraic special frequencies of the black holes?, (the algebraic special frequencies are 
purely imaginary [33] )■ 
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A Why m is a half-integer? 

In Section |21 we use the diagonal triad (j2J). As is well known, it is also possible to employ 
the Cartesian triad [33] for which the Gamma matrices are equal to 

7t (c) = Pli, 

li c) = j (cos(% 2 + sin(0) 73 ) , (96) 
7e C) = r{- sin((9)7 2 + cos(6 l )7 3 ) , 

when the background is 3D de Sitter. In Eq. (J3J), we previously define the matrices 71, 
72, and 73 which appear in Eq. (|96|). The relation between the solutions of the Dirac 
equation calculated in both triads is 03] 

SI = 1(c), (97) 



where 

^ ~ 1 -/simfl/2) c-os(0/2) J ' { 2 
The solution \I/( C ) is single valued [33] , but we note that the matrix S satisfies 

S(6 + 2n) = -S{6), (99) 

therefore the spinor \1/ must satisfy the relation ^(8 + 2tt) = — Since in Section |2] 
the ^-dependence of the spinor \l/ is of the form e tme , we find that the quantity m takes 
the values ±i ±§, ±§, ... 12 



12 Thc author gratefully acknowledges the many helpful suggestions of Dr. G. F. Torres del Castillo on 
the content of this Appendix. 
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B Quasinormal modes in 4D de Sitter spacetime 

In four dimensions, when the background spacetime is spherically symmetric, it is possible 
to show that the equations of motion for a scalar field, an electromagnetic field and an 
odd parity gravitational perturbation simplify to a Schrodinger type differential equation 
[32], For AD de Sitter metric this equation is equal to 



d 2 $ 
dr? 



V{r)]$ = 0, 



(100) 



where r* is the tortoise coordinate defined in Eq. (J34j) and 

+ 2a 



V(r) = P 2 



I 2 



(101) 



with a = 1 for a massless minimally coupled scalar field and a = for an electromag- 
netic field and an odd parity gravitational perturbation. For an even parity gravitational 
perturbation the effective potential V(r) is a more complicated expression 
In de Sitter spacetime Eq. ()100|) becomes 



13 



2\2 



-2z(l-z 2 )- + u 2 -(l-z 2 ) 
dz z dz 



- 2a 



$ = 0, 



where zl = r and u = ul as previously. Making the change of variable y 
dU, [HH we find that Eq. (ITTT21 takes the form 



(102) 

us, na, 



dy- 



+ 



1 



1-2/ 



^ 2 



j(j + l) + 2« 



dy 



%(1 - y) 



cu 2 



j(j + 1) 



$ = 0. 



(103) 



4(1 -y) 2 

Proposing a solution of Eq. ()103|) in the form 

$ = y A (l - y) B $, (104) 
we get that the function $ satisfies a hypergeometric type differential equation 



d 2 d 
y(l - y)— + [c - (a + b + l)y}— - ab } $ = 0, 



(105) 



where the parameters a, b, and c are equal to 



A + B + - + -[l + 8a]2. 

4 4 L 1 ' 

A + B + ---[l + 8a]2. 
4 4 L J 

2A+ -, 



(106) 



13 Nevertheless see Ref. |47j . 
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provided that 



i 4- i 

A = { 2 _J 2 ' B = ±f . (107) 
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In the following we study in detail the case A = ~ + | and B = Thus the solutions 
of the differential equation ()105|) are [2H] 

<§« = 2 F 1 (a,6;c; 2/ ), (108) 
$( 2 ) = ^(a-c+M-c + l^-c;!/), 

because the parameter c is a half-integer. 

Taking $( 2 ) solution of Eq. fjl05jl . we find that the function $ is divergent at the 
origin (y = 0) for j > 1. In the present paper we are studying regular solutions at r = 0, 
hence in the rest of this appendix we only consider Thus 

ill 

$ = y2 + 2(l- 2/ )2 2 F 1 (a,6;c;y). (109) 
If c — a — b is not a integer, we can write Eq. ()1U9|) in the form j2Hj 

$ = ^ ( r^V^ 1 " 2^i(a, b; a + b - c+ 1; 1 - y) 
+ r(a) r( 6 ) ^ ~ y)" T ^i( c - a, c - 6; c - a - 6 + 1; 1 - y) j . (110) 

Since the field $ depends on time as e~* w , it follows that the first term in curly braces in 
Eq. (jllUJl represents an ingoing wave as r — > I and the second term represent an outgoing 
wave as r — > /. As we are calculating the de Sitter QNM frequencies, the field must be 
purely outgoing as r — > /, that is, we must impose the condition 

c — a = — rii, or c—b — —ni. (m) 

Nevertheless, from the expressions for the parameters a, 6, and c given in Eq. (jl06j) . 
it is possible to show that if c — a — b is not an integer, then c — a and c — b cannot be 
integers; therefore we cannot satisfy the conditions (jlllj) . 14 

When c — a — b is a negative integer, 15 that is, c — a — 6 = — n, n = 1, 2, 3, . . . we can 
use the property (|68p of the hypergeometric function to write the field $ as 



i + i £ rr(a4-&-n)r(n) /1 , 

U K y) \ T(a)T(b) 1 U> 



x y- ( a -n) P (b-n) p _ _ (-l) w r(a + b - n) (a) p (b) p _ y 
p\(l-n) p r{a-n)T{b-n) ^p!(n + p)! 

x [In(l - y) - ip(p + 1) - ij}{p + n + 1) + ^(a + p) + ijj{b +p)] } . (112) 



14 We note that (1 4 8a) 2 is an integer for a = and a = 1. 

15 If c — a — b is zero or a positive integer, then there are no QNM frequencies (see footnote |SJ . 
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The first term in curly braces of Eq. (J112)) represents an outgoing wave and the second 
term represents an ingoing wave; therefore the field $ is purely outgoing if we impose the 
condition 

a — n——ni, or b — n = — n\. (H3) 

From Eqs. (jll3|) we find that the QNM frequencies are equal to 

3 1, i 3 \, i . JS 

i£o = j + - + -(l + 8a)2 +2m, iuj = j + - - -(1 + 8a)2 + 2n v (114) 

For an electromagnetic field and an axial gravitational perturbation these frequencies are 
equal to those calculated in Section 13 .2\ since for each j the frequencies ()114j) are even or 
odd integers, depending on the j value. For a massless minimally coupled scalar field our 
result is different from that of Refs. [TB] - 

In our analysis we do not include the mode j = of the massless scalar field. As 
shown in Ref. ^2], when this field is moving in de Sitter spacetime its modes j = and 
j > behave in a different way. 

We can easily write the effective potential given in Eq. (jlOlj) in the form 

V{r) = P {sinh 2 (rj/) " cosh 2 (r»/Z) } ' (U5) 

where r = Ztanh(r*/Z) (r* is the tortoise coordinate). (See Refs. [T3] and [H] for the 
effective potentials of a massive scalar field and a gravitational perturbation, respectively, 
when these fields are propagating in D-dimensional de Sitter spacetime.) 
In Ref. |IH] Beyer proved that for the Poschl- Teller potential 

Vpt(x) = ^° . , x G K, (116) 
cosh (x) 

"after a large enough time to, the solutions of the wave equation corresponding to C°° data 
with compact support can be expanded uniformly in time with respect to the quasinormal 
modes, thereby leading to absolutely convergent series" For r* > 0, is there a similar 
result for the potential (jll5j) ? We believe that this point must be studied further. 



C Otchik's method 



In Ref. j2T] Otchik proved that in order to find exact solutions of the massive Dirac 
equation in AD de Sitter metric, it is necessary to solve the system of differential equations 



[l-z 2 )— + iu 
dz 



'1 - z 2 )— -iu> 



dz 



R 2 = (l-z 2 ) 1 i^-iM) /?,. 
R 1 = (l-z 2 )^(^ + iM) n,. 



;ii7) 



where zl = r, uj = ul, M = Ml as above and k = j + |. This system of differential 
equations is equivalent to that given in Eq. (fTSj) (see Appendix of Ref. [17J). 
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To solve the system of differential equations (|117|) . Otchik proposes that the functions 
Ri and R 2 take the form 

R x = i[b 2 (l - £)( K+1 )/ 2 5( 1 -^)/ 2 / 1 + bx(l - xY /2 x' iG),2 f 2 ], 

R 2 = 6i(l -x)^ 1 )/ 2 ^ 1 -^/ 2 /^ 62(1 -x)^"* 572 ^, (US) 

where x = 1 — z 2 , and the functions b\ and b 2 are equal to 

~ bl = [x 1 /2 + ,(l_ £ )l/2 ] l/2 + ^~l/2_ z(1 _ 5)l/2]l/ 2j 

6 2 = _[£l/2 +i (l_ £ )l/3]l/2 + i[£l/2_ (119) 

Substituting Eqs. (|118|) and (|119|) into Eqs. (|117|) . we find that the functions f\ and 
/2 satisfy the coupled differential equations 

5(1 - xj^± - - ~)(1 - x) + {k + \)x\h = i(- K -\ + i£o + iM)f 2 , 



2 



da; 



w-i-iw + iM)/!. (120) 



From Eqs. fjl2()jl . after eliminating f\, it is possible to show that the function f 2 must be 
a solution of the differential equation 

x(l - x) ( ^rJr + + |) - (k — iuj + l)x]^ 
ax 2 ax 

- \{k -iu-iM + l){K-iu + iM - i)/ a = 0, (121) 
which is a hypergeometric type differential equation. 
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